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Abstract 

A general explicit coupling for mutual synchronization of two arbitrary identical continuous 
systems is proposed. The synchronization is proved analytically. The coupling is given for all 19 
systems from Sprott's collection. For one of the systems the numerical results are shown in detail. 
The method could be adopted for the teaching of the topic. 
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1. INTRODUCTION 



Synchronization is a fascinating phenomenon in nature 0, 0] and a useful one in human 
activity 0]. Very simple systems like driven gravitational pendulums or driven nonlinear 
LCR electric circuits can have chaotic behavior. This means that such systems can not be 
kept to oscillate in synchrony due to their sensitivity to initial conditions. To synchronize 
them an additive coupling have been tried that is feasible from the engineering point of 
view: a constant multiplication of the differences of the states .This type of driving does 
not work in general. Highly elaborated and mathematically based methods are needed 0, 
||| 0. If the above mentioned simple chaotic systems can be synchronized with a precise 
coupling that is not intuitive at all then the coupling of biological cells that are multivariables 
nonlinear systems (like millions of neurons that fire together to control our breathing or 
the coordinated firing of thsousands of pacemakers cells in our hearts yj) is hard to be 
imagined. There are known many results on several types of synchronization Mutual 
synchronization is implied in the synchronization of networks Here we give, for the first 

time, a precise general coupling between two arbitrary identical oscillators in order to get 
synchronization. The proposed coupling is mathematically based and the synchronization 
is analytically proved. Numerical results are shown. The paper is organized as follows. In 
section |21 a coupling is proposed and a theorem is proved for the synchronization of the 
two systems. Section |3] contains detailed calculations for the Hurwitz matrix and numerical 
results for one of the systems. The paper ends with conclusions. 



2. MUTUAL SYNCHRONIZATION 

Let's consider two identical nonautonomous systems of the form: 

x\ = F(x 1 ) + Bu u yi = Cx 1 , x x (0) = x 10 (1) 

x 2 = F(x 2 ) + Bu 2 , y 2 = Cx 2 , x 2 (0) = x 20 (2) 

where Xj G 9ft n , Uj G 3ft n , j = 1,2. For the prescribed inputs u\ = —7(2/1 — 2/2) and 
u 2 = — , y(y 2 — 2/1), the mutual synchronization of (JTJ) and (J2J) has been studied in |9| using 
the powerful concept of passivity l3|- Namely, it is proved that there is a constant 7 such 



that for 7 > 7 the systems (JIJ and (J2J) synchronize and their dynamics is bounded. Here 
we consider a particular case with B = C = I (identity matrix). The inputs u% and u 2 
are not prescribed but they result from the condition of synchronization. In order to avoid 
the indices with two figures we adopt the notation: x\ = x, x 2 = y, U\ = S(t)u x (x,y) and 
u 2 = S(t)u y (x, y). In the following we study the coupled systems: 

x = F(x) + S(t)u x (x, y), x(0) = x (3) 

y = F(y)+S(t)u y (x,y), y(0) = y (4) 
where S(t) is or 1 as a switch. The above systems will synchronize if 

lim \\x(t) — y 

t^oo 
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for any x(0),y(0) G © C 3?". The proposed couplings are: 

u y( x ,y) = [ H — ) — 2 — ^ 

where s = 2±S anc l j s an arbitrary constant Hurwitz matrix (a matrix with negative real 
part eigenvalues). With equations and (jUJ) we announce the following Theorem. 

Theorem. With S(t) = 1, the systems J3J) and with the coupling ([3J) and (GJ) will 
synchronize for any x(0) and y(0) with \\x(0) — y(0) \\ small enough. 

Proof. Along with s = we make the notation r = It follows x = s + r and y = s — r. 
Substracting eq. (JU) from eq. Q one has: 

._ F(s + r)-F(s-r) | / dF(s)\ (?) 

We use the Taylor expansions: 

F ( s+r ) = F ( s ) + ^ r + __U r(rj + ... (8) 

f(s . r ,. fW + ^ H + ^ + ... (9) 

Using the first three terms in eq. © and eq. 0, eq. (0) becomes: 



r = Hr (10) 

With i/ a Hurwitz matrix, eq. (jlO)) assures that r(£) — >• for any r(0) for which the 
Taylor expansions eq. (jHJ) and eq. (jHJ) are valid. This means that x(t) — > ?/(t) for any r(0) = 
x(0) — y(0) with ||x(0) — y(0)\\ small enough. If F(x) is quadratic then the Taylor expansions 
have only three terms and in this case the eq (jl(J|) is not anymore an approximation. 

□ 



Here we note a significant difference between master-slave synchronization and mutual 
synchronization. For master-slave synchronization 0, , the error dynamics is described by 
an approximate equation like eq (J1U|) for any nonlinear F(x). For mutual synchronization 
of systems with F(x) a polynomial up to quadratic terms, eq (fTT)|) is exact. This could lead 
to the conclusion that mutual synchronization can be achieved in an easier manner. This 
is not generally true. A deeper analysis is necessary. In spite of the fact that the eq (fTUj) 
assures that the two dynamics will be closer and closer, it is not sure that their dynamics 
are bounded. Even if they are bounded they could have large excursions in the phase space 
that is very undesirable from the engineering point of view. Further studies are needed to 
establish how to choose the parameter p (or pi^) in the matrix H (see Table I) and the 
initial conditions in order to have bounded dynamics of the synchronized systems. Otherwise 
we can manage this by switching on/off the coupling like this: S(t) = 1 when \\x — y\\ < 5 
and 1 1 a; 1 1 < R where R is the radius of a sphere that contains the attractor of the system 
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x = F(x). In Table II (next section), we give the numerical values of the parameter and the 
initial conditions for which the synchronization was verified numerically. 

Matrix H can be chosen in such a manner that the terms from eq. (0) and eq. © to be as 
simple as possible (HEf and this depends on the particular form of F(x). If is constant 

then we can choose = and the corresponding term (u x )ik will be zero. The simplest 

coupling in eq. (0) will be when F(x) contains one nonlinear term and this contains one 
variable. In this case the coupling contains one term (see Table I, below, systems F ,H, I, J, 
L, M ,N, P, Q, S). We apply this strategy to all systems from Sprott's collection Table 
I contains the coupled eq. © and eq. © for all 19 systems from Sprott's collection. 



3. NUMERICAL RESULTS 

We give detailed calculations for the system S from the Sprott's collection (see Table I). 
The system in x is: 

X\ = —X\ + 4x 2 

X 2 = Xi + x\ (11) 
X 3 = l + Xi 

The Jacobian is: 



-1 -4 

1 2s 3 | (12) 
1 



-1 


-4 





1 





P 


1 









dF(s) 
ds 

We can choose matrix H as: 



(13) 

v 1 ° °y 

The characteristic equation is: 

A 3 + aiA 2 + a 2 X + a 3 = (14) 

where 

ai = 1, a 2 = 4, a 3 = 4p. (15) 
The Ruth-Hurwitz conditions (for eq. (|14|)) are: 

ai > 0, aia 2 — a 3 > 0, a 3 > 0. (16) 
This gives for the parameter p the condition < p < 1. In this case we have: 



u\ = u 1 = 

u% = {p- 2s 3 )r 3 J = (p - 2s 3 )(-r 3 ) (17) 
ul = u? = 



The two coupled systems follow now as (see Figure 1 for numerical simulation): 
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TABLE I: The coupled systems and (j2J of the Sprott's collection J±i] with Si = ' and 



r . = xi=Vi w here i = 1, 2, 3 and S(t) = 1. 


System 


system x 


system y 


parameter 


A 


±1 = X 2 

±2 = —XI + X2X3 + ( — 1 — S 3 )r2 + (1 - 52)^3 
£3 = 1 — x\ + pri + 2S2T2 


J/i = V2 

m = -y\ + j/2 j/3 + (-1 - s3){-n) + (1 - s2)(-r\3) 
2/3 = 1 — 2/1 + P(- r 'i) + 2s2(— n) 


-1 < p < 


B 


±\ = x 2 x 3 + (1 - s 3 )r 2 + (pi - s 2 )r 3 
±2 = XI — X2 

±3 = 1- X\X2 + (P2 + S2)ri + (1 + Sl)r2 


yi = I/2I/3 + (1 - s 3 )(-r 2 ) + (pi - s 2 )(-r 3 ) 
2/2 = 2/1- j/2 

2/3 = 1 - 2/12/2 + (P2 + S2)(-ri) + (1 + Sl)(-T2) 


Pi > 1 
P2 < -1 


C 


± 1 = X2X3 + (p - s 3 )r 2 - (1 + s 2 )r 3 

±2 = XI — X2 

±3 = l-xj + (l + 2si)n 


2/1 = 2/22/3 + (P - S3)(-T2) - (1 + s 2 )(-r 3 ) 
2/2 = 2/1- 2/2 

2/3 = l-2/i +(l + 2si)(-n) 


p < 


D 


±\ = —X2 
±2 = Xi + X3 

±3 = X1X3 + Zx\ + (p2 - sz)r\— 
-6s2»*2 + (pi - si)r 3 


2/1 = -2/2 
2/2 = 2/1 + 2/3 

2/3 = 2/12/3 + + (p 2 - s 3 )(-ri)- 
-6s 2 (-r 2 ) + (pi - si)(-r 3 ) 


T)i < 
Pi < P2 
P2 < 


E 


±\ = X2X3 + (pi - S3)r2 + (p2 — S2)r3 

±2 = x\ — X2 + (1 — 2si)ri 

±3 = 1 — 4x 1 X2 + (-4 + 4s 2 )ri + 4sir 2 


2/1 = 2/22/3 + (pi - S3)(-»"2) + (P2 - S2)(-r 3 ) 

2/2 = y\ - 2/2 + (1 - 2si)(-n) 

2/3 = 1- 4j/ij/ 2 + (-4 + 4s 2 )(-n) + 4si(-r 2 ) 


Pi < 
P2 > 


F 


±1 = X2 + X3 

±2 = —xi + 0.5x2 

X3 = x\ - X3 + (p — 2si)ri 


2/1=2/2 + 2/3 

2/2 = —2/1 + 0.52/2 

2/3 = 2/j? - 2/3 + (p - 2si)(-ri) 


pe (-2,-0.75) 


G 


±1 = 0.4xi + X3 

X2 = X1X3 - X2 — S3Ti + (p - Sl)r 3 
X3 = -XI + X2 


2/i = 0.4yi + y 3 

2/2 = 2/12/3 - 2/2 - s 3 (-ri) + (p - si)(-r 3 ) 
2/3 = -2/1 + 2/2 


p e (-2.5, -0.6) 


H 


xi = -x 2 + x^ + (p - 2s 3 )r 3 
X2 = xi + 0.5X2 
X3 = xi - X3 


2/1 = -2/2 + 2/3 + (P - 2s 3 )(-r 3 ) 
2/2 = 2/1 + 0.52/2 
2/3 = 2/1- 2/3 


pe (-2,-0.75) 


I 


xi = —0.2x2 

X2 = Xi + X3 

X3 = xi + x| - X3 + (p - 2s2)r2 


2/i = -0.2j/ 2 
2/2 = xi + J/3 

2/3 = xi + 2/| - 2/3 + (p - 2s 2 )(-r 2 ) 


p < -0.2 


J 


xi = 2x 3 

X2 = —2X2 + X3 

±3 = -xi + X2 + x\ + (p — 1 - 2s2)r2 


2/1 = 22/3 

2/2 = -2?/2 + 2/3 

2/3 = -2/1 + 2/2 + 2/| + (p - 1 - 2s 2 )(-r 2 ) 


p < 


K 


Xl = XlX 2 - X 3 - S 2 Tl + (P - Sl)»"2 
X2 = Xi — X2 

X3 = xi + 0.3x3 


2/1 = 2/12/2 - 2/3 - S2(-n) + (P - si)(-r 2 ) 
2/2 = 2/1- 2/2 
2/3 = 2/1 + 0.3?/3 


p e (-3.3, -0.51) 


L 


xi = x 2 + 3.9x 3 

x 2 = 0.9xJ - x 2 + (p - 1.8si)ri 
X3 = 1 - Xi 


2/1=2/2 + 3.9?/3 

2/2 = 0.9yl -y 2 + {p- 1.8«i)(-n) 
2/3 = 1-2/1 


p < 


M 


Xl = -X3 

X2 = — x\ — X2 + (p + 2si)ri 
x 3 = 1.7 + 1.7xi + x 2 


2/1 = -2/3 

2/2 = -2/f - 2/2 + (p + 2si)(-ri) 
2/3 = 1-7+ 1.7jyi +2/2 


p e (-1.7,0) 


N 


xi = —2x2 

X2 = xi + x\ + (p — 2s3)r3 
x 3 = 1 + x 2 — 2x 3 


2/1 = -2?/2 

2/2 = 2/1 + 2/3 + (P - 2s 3 )(-r 3 ) 
2/3 = 1 + 2/2- 2y 3 


p < 


O 


Xl = X2 

X2 = Xi — X3 

X3 = xi + X1X3 + 2.7x'2 - S3n + (p - si)r 3 


2/1 = 2/2 

2/2 = 2/1- 2/3 

2/3 = 2/1 + 2/12/3 + 2.7j/2 - s 3 (-ri) + (p - si)(-r 3 ) 


pe (-1,-0.37) 


P 


xi = 2.7x2 + X3 

X2 = -xi + x\ + (p - 2s2)r2 

x 3 = Xl + x 2 


?)i = 2.7w> + ?/3 

2/2 = -2/1 + 2/| + (p - 2s 2 )(-r 2 ) 
2/3 = 2/1 + 2/2 


P€(-I,i) 


Q 


xi = -x 3 
X2 = Xl - x 2 

x 3 = 3.1xi + x\ + 0.5x 3 + (p - 2s 2 )r 2 


2/1 = -2/3 
2/2 = 2/1" 2/2 

2/3 = 3+2/1 + y\ + 0.52/s + (p - 2s 2 )(-r 2 ) 


pe (-3.1,-1.8) 


R 


xi = 0.9 — X2 
X2 = 0.4 + X3 

x 3 = xix 2 - x 3 + (1 - s 2 )ri + (p - si)r 2 


2/1 = 0.9 - J/2 
2/2 = 0.4 + j/3 

2/3 = 2/12/2 - 2/3 + (1 - s 2 ){-n) + (p - si)(-r2) 


p < -1 


S 


xi = — xi — 4x2 

X2 = xi + x\ + (p — 2s3)r3 

x 3 = 1 + X'l 


2/1 = -2/1 - 4 2/2 

2/2 = 2/1 + 2/3 + (P - 2s 3 )(-r 3 ) 

2/3 = 1 + 2/1 


pe (0,1) 
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FIG. 1: Numerical results. S(t) = for t < 120 and S(t) = 1 for t > 120 for the coupled systems 
((THJ) and (JTJ3). Initial conditions are zi(0) = x 2 (0) = x 3 (0) = -1 and yi(0) = y 2 (0) = y 3 (0) = 1. 
Top figure represents x% and y\ and the bottom figure is the error function x\ —y\. It takes about 
15 time units since the mutual coupling is switched on (at t = 120) until the synchronization is 
achieved. 



x\ = —X\ + 4^2 

x 2 = Xl + 4x 2 3 + - 2s 3 )r 3 (18) 

X 3 = 1 + Xi 



y\ = -yi + 4y 2 

2/2 = 2/1+ 4i/3 2 + S(t)(p - 2s 3 )(-r 3 ) (19) 
Va = l + 2/i 



where S3 = X3 + V3 and r 3 = X3 ~ m . 

In the same manner can be written the couplings for all 19 systems from Sprott's collec- 
tion. If the conditions ()16|) can not be fulfilled with one parameter p then a second parameter 
should be introduced (see Table I, systems B, C, D, E, O, R). The Hurwitz matrix should 
be chosen the same as in 0. The choice of H is not unique for systems that have several 
terms in the coupling term. And, as was mentioned in the previous section, the theorem 
does not assures that the dynamics of the coupled systems is bounded. 

In Table II we present numerical values for the parameter p or pi,p 2 (see Table I) and 
initial conditions for which the synchronization was verified numerically. It can be observed 
that it was not difficult to find such numerical values. They are rather homogeneous. When 



6 



TABLE II: Numerical values of the parameter p or p±,P2 and initial conditions for the systems in 
Table I for which the synchronization was verified numerically. 



System 


A 


B 


C 


D 


E 


F 


G 


H 


I 


J 


K 


L 


M 


N 


O 


P 


Q 


R 


S 


Xi(0) = - yi (0) 
i = 1,2,3 


1 


1 


1 


0.1 


1 


0.1 


0.1 


1 


1 


1 


1 


1 


0.1 


1 


0.1 


0.1 


1 


1 


1 


P 


-0.1 




-1 






-1 


-1 


-1 


-10 


-1 


-1 


-1 


-1 


-1 


-0.5 


-0.5 


-2 


-2 


0.5 


Pi 
P2 




2 

-2 




-0.1 
-0.05 


-1 
1 































the dynamics was far from the dynamics of the original system x = F(x) for x(0) = — y(0) = 
(1,1,1) then we changed to the initial conditions x(0) = — 2/(0) = (0.1,0.1,0.1). Further 
studies are needed to find analytic conditions for the matrix Hurwitz H and the initial 
conditions that assure a bounded dynamics for the coupled system. 

Using the same notation like the Table I and the Hurwitz matrix found in jfj, two Lorenz 
systems mutually synchronized look like: 

±i = s(-x l + x 2 ) 

x 2 = rxi - x 2 - xix 3 + S(i)(sxX 3 + (p + s 3 )ri) (20) 
x 3 = -bx 3 + x x x 2 - S(t)(s 1 r 2 + s 2 ri) 

Vi = s(~Vi + 2/2) 

2/2 = ryi -y 2 - ym + S(t)(-siys + (p + s 3 )(-ri)) (21) 
m = ~by 3 + ym + S(t)(sir 2 + s 2 ri) 

For (s,r,b,p) = (16, 45.6, 4, -60) and x(0) = — 2/(0) = (1, 1, 1), the synchronization was 
verified numerically. In Q, the mutual synchronization of two Lorenz systems was obtained 
by using one term in the first equations and any intial conditions. 

Also for the simplest chaotic sytems x + ax — x 2 + x = 0, the coupled systems that will 
synchronize look like: x +ax — x 2 + x+ (p+ (x+y))(x — y) = and V +ay — y 2 + y + (p+(x + 
y))(-x+y) = and for (a, p) = (2, -1) and (x(0), i(0), x(0)) = -(y(0), y(0), y(0)) = (1, 1, 1) 
the synchronization was checked numerically. 

4. CONCLUSIONS 

A precise analytical scheme for mutual synchronization is proposed. The method is 
mathematical rigorous in terms of choosing the expression of coupling and indications are 
given how to avoid unbounded dynamics of the coupled systems. 
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